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We discuss the role of striction in the intertwined magnetic and structural phase transitions in 
the underdoped iron-pnictides. The magneto-elastic coupling to acoustic modes is then derived and 
estimated in framework of the multiband spectrum for itinerant electrons with nesting features. We 
argue that the 1-st order character of the magneto-elastic phase transition originates from the lattice 
instabilities near the onset of spin-density wave order introducing, thus, a shear acoustic mode as a 
new order parameter. Taking non-harmonic termis in the lattice energy into account may explain 
the splitting of the structural and magnetic transitions in some oxypnictides. Fluctuations of the 
magnetic order parameter show up in the precursory temperature dependence of the elastic moduli. 



PACS numbers: 74.70.-b,74.20.-z,74.20.Rp,74.20.Mn 

Considerable efforts have been concentrated recently 
on studies of the interdependence for the antiferromag- 
netic (AFM) and the superconducting (SC) transitions 
in the new iron-based compounds Doping is known 
to destroy AFM and increase Tc - the temperature for 
onset of SC0. 

In this presentation we address the pecuharities and 
interrelations for the magnetic (Tm) and structural 
(Tstr) transitions in the parent and underdoped iron- 
pnictides. For concreteness, we discuss below the prop- 
erties of two systems only: the quarternary, REFeAsO 
("1111"; RE stands for a rare earth), and the bi-layered, 
AFe2As2("A22", where A=Sr,Ca,Ba), materials. In both 
classes the intertwinned magnetic and the structural 
transitions of the weak 1-st order are observed in a tem- 
perature interval of ^ 100-200 K. In llll's the structural 
change precedes the magnetic ordering 0, while in 
A22's magnetic order and change in the lattice symmetr 
occur simultaneously in the single 1-st order transition 

It is known that coupling to the lattice may trans- 
form a magnetically driven transition into the weak 1-st 
order transition accompanied by structural changes jsj. 
In Ref. [6,] the problem has been rigorously solved for 
the elastically isotropic solid, taking fluctuations into ac- 
count. Unfortunately, the method of Ref. 6] does not ap- 
ply to anisotropic materials, such as iron-pnictides with 
the tetragonal symmetry and layered structure. The so- 
lution of Ref. however, undoubtedly contains the main 
physics for these phenomena. Correspondingly, below we 
simplify the approach of Ref. @ to model the magneto- 
elastic interactions (striction) and phase transitions in 
the parent and underdoped FcAs-systems. 

The magnetic order in oxypnictides is built of the al- 
ternating (AFM) spins running along one of the lattice 
axis, while in the other direction the spins are ordered 
ferromagneticallyQ. In the tetragonal notations, the 
structure vectors are commensurate: Qi — (0, tt) and 
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Q2 = (tt, 0). With the local picture of interacting Fe 
spins and the two Heisenberg exchange interactions con- 
stants Ji (nn) and J2 (nnn), such a ground state realizes 
itself at the inequality J2 > Ji/2. Accompanying orthog- 
onal lattice distortions can be understood, as in Ref. [5|, 
in terms of a "spin-Peierls" effect, i.e., the variation of 
the exchange integrals at the lattice deformation [?]. 

The new materials, however, are semimetals. Pnic- 
tides are better described in terms of an itinerant scheme. 
The consensus is that the energy spectrum obtained 
in the "first principle" calculations [8] presents a good 
starting point. According to 8], the electron spectrum 
bears the multiband character: there are Fermi surfaces 
(FS) for the two hole-pockets (h-) at the F-point (0, 0) 
and two electronic pockets (e-) located in the tetragonal 
(unfolded) [il reciprocal lattice at (0, tf) and (tt, 0). An 
approximate nesting between the e- and h-pocket is be- 
lieved to be responsible for a spin-density wave (SDW) 
instability with the two vectors, Qi and Q2, mentioned 
above [9J- The main features of this spectrum have been 
reproduced in many numerical calculations for all classes 
of the new Fe-based materials. In Ref. fio'l this spectrum 
has been directly observed for LaFePO in the de Haas- 
van Alphen (dHvA) experiments. Below we accept this 
model. With ~ 50K, - T^tr ^ 100-200ii: and the 
Fermi energy Ep ~ 0.1 — Q.2eV for the pockets' sizes, the 
model is expected to allow a mean field treatment both 
for SC and magnetic phenomena [Tl|. 

"First principle" calculations provide a reasonable de- 
scription of the ground state properties, albeit some dis- 
agreements are not uncommon in the literature (see in 
Ref. [l2|)- However, subtleties, especially in a vicinity of 
phase transitions, remain beyond the reach of numerical 
analysis. 

Among advantages of the "nesting" modeljl, M, H 
is that its formalism is practically identical to the 
well-studied BCS scheme for SC. It is a weak-coupling 
mean- field scheme if Tsdw ^ Ep- (Fluctuations be- 
come important only in a narrow temperature interval, 
I AT I /To = G« ^ 1, as given by sort of a Levanyuk- 
Ginzburg criteria, Gi.) The Landau functional near the 
transition can be derived exactly as it has been done for 
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A(p) = 



A(p) + (z^(p) terms ) 



FIG. 1: (a) The logarithmic nesting model (b) Mean field 
equations near Tc- 



superconductivity near T^|l£ 



Nevertheless we prefer not to write equations explicitly. 
Numerous unknown parameters that include parameters 
of the e- and h-pockets, interactions and six elastic mod- 
uli for the tetragonal symmetry of pnictides should make 
it useless. For the sake of transparency we hold, where 
possible, the discussion on the qualitative level. 



The mathematical analogy of the model (Tsl . [l4| with 
the Cooper pairing in a BCS-like scheme stems from 
a logarithmic divergence for the scattering in the e-h- 
channel in Fig. [T^, which has the following form: 
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= Vi^{EF)\n {Cj/T,6), 



(1) 
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where 



(2) 



are e- and h- energies, while S accounts for deviations 
from the ideal nesting. With the notations, Ai and A2 
for the triplet (SDW) order parameters corresponding 
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to the structure vectors Qi = (0;7r) and Q2 = (tt, 0), 
the equations for A-s have the same structure as in SC 
theory [ill, as demonstrated in Fig. [T)d. Let Tq be the 
"bear" transition temperature in the absence of coupling 
to the lattice. The Landau functional near Tq is of 
the form: 
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(3) 



where v = v{Ef) is the characteristic density of states 
(DOS). A, cq, 60 a-re constants, and t ~ [T — Tq)/Tq; 
Co — vf/^ttTq. The functional Eq.® is the starting 
point in the general theory of the second order phase 
transitions [TgI. \it\ . To account for interactions with the 
lattice, one adds to Eq.((3]) a coupling to the lattice of the 
form: 



(4) 



Hstr = -q / uir)A{rydr, 



and the elastic energy itself that we also write schemati- 
cally as: 



(5) 



where K is an elastic modulus (such as the bulk modulus, 
for instance); the notation u stands for components of the 



strain tensor. 
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duj 
dxk 



dxi ) 



(6) 



For magnetic phenomena, Eq.Q describes magneto- 
elastic interactions and q is called the striction constant. 
The total Gibbs energy near the phase transition is 



He 



H, 



(7) 



In the itinerant model of Eq.((T]) the parameters Ai and 
A2 are built on the Bloch wavefunctions. These spa- 
tial features present in the SDW parameters immediately 
lead to the quadratic coupling between A-s and acoustic 
degrees of freedom {Q = 0) given by Eq.(|4]). An esti- 
mate of the striction constants in Eq.(I3]) follows directly 
from the logarithmic contribution of Eq. ([T|) . Indeed, elas- 
tic deformations change the parameter 5 in Eq.® that 
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controls the degree of the nesting: 

S = S{u) = So + Xu 
From Eq.([T|) and Fig. [TJa it is seen that 

u 



A- 



Tn 



Making use of Eq.Q, we estimate 
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(8) 



(9) 



(10) 



(taking A of order of characteristic atomic energy, Ep). 

Return to Eqs. (H))-© and assume homogeneous u. 
Minimizing $ over u, we obtain: 
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(11) 



(with V for the volume). The effective Landau functional 
(l is 



1 

2K 



jA^ir)dr 



V 



(12) 



In the mean field A^(r) = const, and Eq. p^ repro- 
duces the result of Ref.Q: the functional, fi, acquires a 
negative contribution to the biquadratic term in Eq.(I3]). 
Should the total be negative, the transition becomes 1- 
st order. (Higher order terms in the expansion Eq.([3]) 
become necessary.) 

In the presence of an external homogeneous deforma- 
tion, iiext substitution into Eq.®, r — > t + X{uext/To), 
leads to: 
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^ext 



Differentiating, we find: 



A" d^n 



To 



(13) 



(14) 



where C(t) is the specific heat. According to Eq. (fTi)l . al- 
ready in the Landau mean field approach the jump AC 
in the specific heat at the transition is accompanied by 
a negative step-like change of elastic moduli. The varia- 
tion is of the atomic order. Therefore the lattice would 
become unstable, and the 2-nd order character of a tran- 
sition changes to the 1-st order one. 

In Refs.jy, [l^, the attention has been drawn to the 
fact that the specific heat C(T) in Eq. ([14]) actually has 
a singularity at Tq: 



CiT) 



(15) 



[a is the so-called scaling index for the specific heat, 
C{T).) As to the order parameter, it appears below Tq 
(see Ref.dl, §148, p. 484): 



A(T) oc {-Tf{T < 0) 



(16) 



Therefore an elastic instability is ubiquitous for any 
magnetic transition. In the model chosen above the in- 
stability occurs inside a narrow temperature interval con- 
trolled by the Levanyuk-Ginzburg parameter. Provided 
that layered pnictides can be treated, in the first approx- 
imation, as two-dimensional (2D), we write: 



AT 



To 



= G: 



2D 



f TsDW 



< 1. 



(17) 



By the order of magnitude, the same parameter also con- 
trols the width of the hysteresis. 

Note that even for homogeneous lattice deformations 
the magneto-elastic coupling, Eq.(|l]), leads to Eg. p^ . 
with the new biquadratic term that describes non-local 
interactions for the order parameter A. In Ref. inho- 
mogeneous fluctuations result in the cancellation of non- 
local terms Eq. (fT2|l for the isotropic (liquid) media. How- 
ever, the non-local coupling terms are always present for 
any anisotropic solid. As it was mentioned above, the 
problem could be rigorously solved in Ref[fi] only for the 
elastically isotropic solid (i.e., for a solid characterized by 
the bulk and shear moduli). In all other cases the non- 
local terms Eq. that come about due to exchange by 
acoustic phonons become strongly anisotropic. We now 
simplify the model Q by restricting our consideration to 
homogeneous lattice fluctuations. Inhomogeneous fluc- 
tuations should not qualitatively change the physics of 
the problem. 

With this in mind, one can rewrite Eq.([3]) in the form: 



K , 
2 



n 



X 



u — au, 



(18) 



where + ^ is the exact functional Eq.®, i.e., 

the functional that would describe the second order phase 
transition driven by the order parameter A neglecting 
striction effects. In particular, it has the contribution, 
^smg: responsible for the singular behavior of the speciflc 
heat, C(T), Eq. (fT5|) . The term —aii stands for applied 
external stress (a is a proper component of the stress ten- 
sor; this term is needed at calculations of the temperature 
behavior of elastic moduli across the transition.) 

In the presence of non-zero stress, (T ^ 0, fluctuations 
take place around the new equilibrium point, iiext- 



'^ext 



K 



Rewriting in Eq. (jl8p u ^ u + iiext i we obtain 
^2 sr ^ ^ X& 



(19) 



(20) 



The contribution Xa / {TqK) in the argument of ^{x), 
Eq. ([20]) . shifts the transition temperature To under ap- 
plied stress. 



To (a) = To 



Act 

'k 



(21) 
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The Gibbs energy Ea. ((20|) must be minimized over fluc- 
tuations of (homogeneous) deformations, u: 



(22) 



where const is a numeric factor that depends on the 
model details. (Strictly speaking, the values of B in 
Eq.(l27l) differ by factor 2^/2 on the two sides of 
(see §146, footnote on p. 475).) For the singularity in 
the specific heat. Eg. ([27]) gives: 



where 



To 



\a_ 



(23) 



One may easily recognize in Eqs ^U\\ and the anal- 
ogy to equations in Ref. Q describing the transitions in 
the elastically isotropic solid body. 

To be now more specific, note that the elastic energy 
of homogeneous deformations for the tetragonal lattice 
can be rewritten in terms of the even irreducible repre- 
sentation of the Difi group: 

-^Ig • '^xx'^'^yy : ^zzj B\g . Uxx ^yy: ^2g • ^xy^ • ^zx^ '^zy- 

(24) 

In principle, any distortion Eq. (l24|) perturbs the elec- 
tronic spectrum and, hence, affects the nesting features. 
Therefore, the number of the independent striction con- 
stants, g-s, in Eq.(|l]) (or A-s in Eqs ([5]), (O) may be large. 
In view that the layered character of pnictides makes 
them close to two-dimensional systems, the nesting pa- 
rameter S in Eq.® is controlled mainly by the strain 
components in Eg. ipi)) that do not have the z-indicies. 
Rewriting the 2D-part of the elastic energy Eq. ^ in the 
form 



K 



{Ux 



Ml 



'•yy! 



{Ux 



'•yy) 



2m2< (25) 



5C{t) = —B\t\ 

4Jo 



-1/2 



(28) 



Comparing Eq. ((28|) with the normal specific heat, C„ 
vTo , leads to the criterion: 



B 



(29) 



(In the strictly 2D limit (^^o ^ oo) one would obtain 
the criteria Eq.!!?])). With Eqs HI]) and one finds 
that corrections to the elastic moduli of Eq. (P5|) become 
strong in the same order of magnitude temperature in- 
terval Eq.([T7l) (or Eq.(l29l)) if both K and ^ are of the 
atomic scale. 

To demonstrate the emergence of the 1-st order tran- 
sition from Eqs (pni) . (|22p . assume, for simplicity sake, 
a = 1/2 in Eg. lfTS)) . i.e., extrapolate Eg. ((27)l over the 
whole fluctuation interval. From Eq. ([22|) (for each u+,-) 
it follows: 



2 To \x\ 2 To \x\ 



With the notations 



a4 



3 A2 3 A^ 



2-T2' 



2~T§ 



1/2^ 

(30) 
(31) 



with three independent moduli according to Eq. (|24p and 
three magneto-elastic constants qi (or A^ in Eqs. ([8]), ^) 
and taking into account the uniaxial symmetry along 
(0,7r) or (tt, 0) directions for the parameters Ai and A2, 
we simplify the problem further and rewrite Eq.Q as 

T T + ^ {Uxx + Uyy) + ^ (U:r2; - Uyy) (26) 

Jo -to 

The second term is responsible for the orthorombic de- 
formation of the lattice. 

For the actual calculations one would need the expres- 
sion for ^{t). In our case the fluctuations are strong only 
in a narrow vicinity of the phase transition, and we ex- 
pect that the 1-st order transition occurs inside the same 
interval. Therefore only the singular part, 5^lsing{T)^ is 
of the importance. We limit ourselves by the first fluctu- 
ation correction to the mean field O, which we calculate 
in exactly the same manner as in Ref.[T9] (see §147, prob- 
lem, p. 482). The minor difference is that for the strongly 
anisotropic pnictides one needs to introduce in Eq. ^ the 
in-plane and out-of-plane coherence lengths, ^o|| ^ ^O-L- 
After an elementary calculation it follows: 



5Vtsing{T) = -const(rofQ||^4j_i) 



|3/2 _ 



-B|t|3/2, (27) 



\+u+ X+U+ _ 

x = T + p + fj., (32) 

-to -to 
we rewrite Egs. ipOl) as 



a+ \x\ a- \x\ 



or 



a^^p = a_^fi; x = r + p ( 1 H ) = t + v. (34) 



(33) 



Finally, with 



V = {a+ + a-)- — r x^ 



we arrive to the single equation 



X = T + (q!+ -f Q!-)-j — r X 



(35) 



(36) 



Eq. ([36|) like Eq.(13) of Ref. 0, reveals the typical features 
of a 1-st order transition: at small enough r there are 
three solutions for x. If coefficients B in Eg. (ITT)) are 
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equal on both sides of Tq, the transition takes place at 
T = 0, where 



in addition to the elastic terms: 



(37) 



determines the jumps of at the transition ( note 

that at negative X- ^ Q the driving parameter, according 
to Eq. pB]) . is finite). The area of hysteresis is determined 
by the equation: 

^1/2 



2 = (a+ + a^)\x\ 



(38) 



These results qualitatively agree with the observed si- 
multaneous onset both of orthorombic distortions and 
the "stripe" SDW order in the iron system A22yi]. The 
distortions Eg. ([5^ just accompany the magnetic transi- 
tion. In case of the "1111" class the temperature, Tgtr, 
for the onset of the orthorombic deformation precedes 
onset of the magnetic order at Tm [3] . The two tempera- 
tures are rather close: AT = Tstr — Tm > is of order of 
10 - 2QK, so that by order of magnitude AT/Tstr ^ 0.1 
falls into the range of Eq. ^ with Ep ^ O-l - 0.2eV 
taken for the pockets' depths obtained in "first princi- 
ple" calculations d [III . Such closeness seen among most 
of REFeAsO (e.g. compare LaQ and Nd [l^l) is strongly 
in favor that the structural transition is directly related 
to the magnetic instability. In the language of local Fe- 
spins the attempt was made in [2lj to ascribe the tem- 
perature interval separating Tstr and Tm to appearance 
of a "nematic phase" that comes about due to strong 
spin fluctuations above Tm- We suggest that both tran- 
sitions have the common origin and come about as the 
result of the lattice instabilities caused by striction. In- 
deed, according to Eq. fHl) . striction triggers softening of 
elastic moduli as temperature approaches the transition 
interval Eg. p?)) or Eq. (|29|) . Assume that it takes place 
more strongly for modulus m, in Eq.([25l). Recall that 
the orthorombic distortion is the symmetry change by it- 
self and for the tetragonal lattice is characterized by the 
symmetry parameter u_ = Uxx — Uyy So far as the de- 
pendence SflsingiT + {X- /Tq)u^) OH u_ is the only form 
of the nonlinear elastic energy, the above analysis ap- 
plies. However, when the renormalized modulus, jJ-ieff^ 
becomes small, other non-linear terms ever present in the 
lattice must be also taken into account. As the result, the 
Landau functional for the parameter u_ = Uxx ~Uyy also 
depends on those contributions. The higher order terms 
in U- becoming important when renormalized Hieff is 
small nearby magnetic Tq. The mean field treatment of 
the new symmetry parameter, u_, could then be applied 
in the usual way for the 2-nd order tetra-ortho transition 

at Tstr- 

Experimentally, the structural distortions at Tstr ap- 
pear in the weak 1-st order transition|[3|]. That last re- 
sult immediately follows from the fact that by symmetry 
the orthorombic transition in the tetragonal lattice infers 
its own quadratic striction- Indeed, the following cubic 
terms are allowed in the tetragonal lattice by the sym- 
metry reasons: 

(7||7i+u^ + g^UjxU^ (39) 



K 



(40) 



In terms of our model Eqs ([H),©,© the fact that Tm < 
Tstr should mean that the non-zero distortion m_ below 
Tstr makes the nesting conditions Eq.([T]) worse. 

Let's add a few final comments to the analyses above. 
The order of magnitude estimate of the in-plane lattice 
distortions at low enough temperatures follows from Eqs 



K 



vE 



KTn 



F rr.2 

J n 



Ef 



IQ-^ - 10 



-3 



(41) 



that agrees well with the experimental data of Ref . |23. 
23, 24!]. The hysteresis, AT, from Eq.® agrees by the 
order of magnitude with data [J| 

There are no experimental data for the volume change, 
AV/V. The lattice deformation along the c-axis, in prin- 
ciple, could be obtained from Eqs ((39|) and (|40)) together: 



QJ- 2 



(42) 



In Eq. (|42l) both q±^ and K± are expected to be small in a 
layered material. These parameters remain unknown in 
the oxypnictides. 

Eq. (fH)) together with Eq. (P7)l provides another observ- 
able feature: at temperatures above Tq the fluctuation 
corrections to the elastic moduli should behave as 



[T - To) 



-1/2 



while in the 2D limit of isolated planes 
^ (T-To)-i. 



(43) 



(44) 



The data[2J] do not allow to distinguish between the ex- 
ponents in Eqs(|121), 63]) ■ 

To summarize, in the frameworks of multiband elec- 
tronic spectrum with nesting features the theoretical 
scheme is elaborated to treat striction in iron-pnictides. 
Magneto-elastic coupling changes the second order char- 
acter of magnetic transition. The transition becomes of 
the 1-st order. Whether the transition bears strong or 
weak 1-st order character may depend on details. Pro- 
vided that To ^ Ep, the weak 1-st order transition is 
predicted. The model, when applied to the layered FeAs 
systems, leads to estimates of the correct order of magni- 
tude. Discontinuities of all parameters at the transition 
are due to lattice instabilities. Magneto-elastic interac- 
tions may split the magnetic (SDW) transition at Tm and 
the orthorombic deformations at Tstr > Tm- The model 
predicts a noticeable precursory temperature dependence 
above the transition temperature in the elastic moduli. 

The obtained results are in good qualitative agreement 
with peculiarities of the phase diagram of new parent or 
underdoped FeAs materials well above the temperature 
of superconducting transition. 
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